#include "const_nucl.h"

LEC_CONST::LEC_CONST() {
  this->LEC_gA = 1.29;

  int flag = 1;
  /*----------------------------------------------------------------------*/
  if (flag == 1) {
    /*  500 MeV */
    /* PRL 99, 042501 (2007) When CD = - 1.0 CE = - 0.34 From figure FIG. 1 */
    this->LEC_CD = -1.0;
    this->LEC_CE = -0.34;
    this->LEC_C1 = -0.81; //!! GeV^-1
    this->LEC_C3 = -3.2;  //!! GeV^-1
    this->LEC_C4 = 5.4;   //!! GeV^-1//
  }
  /*----------------------------------------------------------------------*/

  /*---------------------------------------------------------------------*/
  if (flag == 2) {
    /* New LECs from Luigi */
    /* 414 MeV
    >> c1=-0.81
    >> c2=3.28
    >> c3=-3.0
    >> c4=3.40
    >> cD=1.85    cD and cE have been fitted on triton binding energy and beta
    >> decay
    >> cE=0.313
    */
    this->LEC_C1 = -0.81; //!! GeV^-1
    this->LEC_C3 = -3.0;  //!! GeV^-1
    this->LEC_C4 = 3.4;   //!! GeV^-1
    this->LEC_CD = 1.85;
    this->LEC_CE = 0.313;
  }
  //
  /*----------------------------------------------------------------------*/
  /*  500 MeV */
  /*  CD = 0.95 CE = -0.041  */
  /* by Laura Marcucci so to obtain the binding energy and beta-decay lifetime */
  /* of triton as the experimental one (or at least very close to them). */
  if (flag == 3) {
    this->LEC_CD = 0.95;
    this->LEC_CE = -0.041;
    this->LEC_C1 = -0.81; //!! GeV^-1
    this->LEC_C3 = -3.2;  //!! GeV^-1
    this->LEC_C4 = 5.4;   //!! GeV^-1//
  }
  /*----------------------------------------------------------------------*/

  /*  450 MeV  Nuclear Matter */
  if (flag == 4) {
    /*  450 MeV  Nuclear Matter */
    /*  power of regulator : 2  */
    /*  https://arxiv.org/abs/1807.06640  */
    this->LEC_C1 = -0.74; //!! GeV^-1
    this->LEC_C3 = -3.61; //!! GeV^-1
    this->LEC_C4 = 2.44;  //!! GeV^-1//
    this->LEC_CD = 0.935;
    this->LEC_CE = 0.12;
  }
  /*----------------------------------------------------------------------*/
  /*----------------------------------------------------------------------*/
  /*  500 MeV  Nuclear Matter */
  if (flag == 5) {
    /*  500 MeV  Nuclear Matter */
    /*  power of regulator : 2  */
    /*  https://arxiv.org/abs/1807.06640  */
    this->LEC_C1 = -0.74; //!! GeV^-1
    this->LEC_C3 = -3.61; //!! GeV^-1
    this->LEC_C4 = 2.44;  //!! GeV^-1//
    this->LEC_CD = 0.495;
    this->LEC_CE = -0.07;
  }
  /*----------------------------------------------------------------------*/
  /*----------------------------------------------------------------------*/
  if (flag == 6) {
    /*https://journals.aps.org/prc/pdf/10.1103/PhysRevC.91.051301*/
    /* NNLO sat  */
    this->LEC_CD = 0.816806;
    this->LEC_CE = -0.039575;
    this->LEC_C1 = -1.121521; //!! GeV^-1
    this->LEC_C3 = -3.925006; //!! GeV^-1
    this->LEC_C4 = 3.765687;  //!! GeV^-1
  }
}
